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We calculate Auger spectra given by the two-hole Green's function from orbitally degenerate 
Hubbard-like models as a function of correlation strength and band filling. The resulting spec- 
tra are qualitatively different from those obtained from fully-filled singly degenerate models due to 
the presence of screening dynamics and multielectron excitations. Application to a real system shows 
remarkable agreement with experimental results leading to reinterpretation of spectral features. 



More than seventy years after its discovery [Q, Auger 
processes continue to evoke strong research interest in the 
general community |^ , since it is potentially a very pow- 
erful probe to investigate the electronic structure of any 
solid. Briefly, the Auger process involves the decay of a 
core-hole by a nonradiative transition involving two elec- 
trons, one filling the core-hole and the second one picking 
up the excess energy and making a transition into the 
continuum. The initial state of the normal Auger tran- 
sition can be considered to be a singly ionized core-hole 
state and the final state is a doubly ionized state together 
with an electron in the continuum. When a core-hole (C) 
decays via the Auger process involving two valence (V) 
electrons, it is usually termed a C-VV process. One can 
anticipate the C-VV Auger spectral shape consider- 
ing two limiting scenarios. The two final state holes may 
be found on different atomic sites due to the presence of 
finite hopping interactions. In such a case they will be 
screened from each other and the Auger process can be 
described essentially by the self-convolution of the occu- 
pied part of the valence band density of states (DOS). 
In the opposite scenario, strong correlation effects may 
localize the holes at the same atomic site, shifting the 
corresponding energy higher than the screened two-hole 
case by approximately the Coulomb interaction strength, 
U . These expectations found a theoretical basis in the 
work of Cini [ || and Sawatzky |^ who calculated the two- 
hole Green's function for fully-filled singly degenerate 
Hubbard-like models. The results indeed showed these 
two spectral signatures, often termed the correlated (or 
the satellite) and the uncorrelated (or the main) features. 
According to such analysis. Auger spectroscopy has the 
unique advantage of providing a direct experimental esti- 
mate of U and such an approach has often been adopted 
for the very same purpose 

The assumptions in the Cini-Sawatzky (CS) theory, 
however, may limit its applicability to real systems. The 
model assumes one orbital per atomic site, whereas most 
real systems invariably involve orbital degeneracies. The 
other assumption of a fully-filled band limits the use of 
the model to few real systems such as Cu and Zn. It is not 
a-priori evident what the consequences of relaxing these 



limitations would be on the Auger spectra compared to 
those inferred from the CS theory. Therefore, we theo- 
retically investigate the Auger spectra from a multiband 
model as a function of band filling and Coulomb interac- 
tion strength in order to understand real systems better 
and to extend the usefulness of this very powerful spec- 
troscopic technique. We find that the results are qual- 
itatively different for partially filled bands compared to 
the fully-filled case. These strong changes occur primar- 
ily due to the screening response to the core-hole in the 
initial state and multielectron excitations in both initial 
and final states. While these processes are not possible 
for a fully-filled case, our results clearly show that such 
processes are in fact the dominant ones in determining 
the Auger spectral shapes in the most often encountered 
partially filled systems. 

The specific multiband model that we consider is based 
on a regular tetrahedron with four atomic sites each con- 
taining triply-degenerate p-orbitals allowing us to treat 
the full multiplet (orbital dependent) electron-electron 
interactions. However, the qualitative features presented 
here are independent of this specific choice, as we have ex- 
plicitly checked by calculating the corresponding spectra 
with three degenerate s orbitals at each site, instead of 
the p-orbitals, as well as with five degenerate d-orbitals. 
The hopping interactions between different orbitals are 
given in terms of Slater-Koster parametrization 0|. In 
order to understand the influence of the core-hole and the 
consequent screening dynamics, we have performed two 
sets of calculations, one with the core-hole potential and 
the other without. The core- hole potential is included 
in the Hamiltonian in the usual way , by including the 
term, -UcTicnd, which effectively lowers the diagonal elec- 
tronic energy by Uc at the core-hole site in the initial state 
of the Auger process. This assumes that the initial state 
is a completely screened core-hole state. This is reason- 
able for the transition elements where the life-time of the 
core-hole is longer than the screening time-scale, but may 
limit the applicability to other systems, such as the lan- 
thanides, where these two time-scales may compete with 
each other. We take Uc to be 1.2 times the effective U 
within the valence band following the usual practice. We 
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present here the results for the number (n) of electrons 
per atomic site being 1, 2, 3, 4, 5 and 6, with the last one 
(n=6) being the fully-filled case. We obtain the ground 
state and the corresponding two-hole Green's function by 
the Lanczos and the modified Lanczos methods 

In Fig. 1 we show some selected two-hole spectra for 
n=6 with a few values of U/W. The Auger spectrum for 
the noninteracting limit {U =0) is the self-convolution of 
the occupied DOS with a width equal to twice the oc- 
cupied part of the single-particle bandwidth, W. With 
increasing U, the intensity within the uncorrelated en- 
ergy region (~ 0-2.2) rapidly decreases, with increasing 
spectral weight thrown out in a narrow feature. This nar- 
row feature with its characteristic dependence on U has 
primarily two-holes at the same site and is easily iden- 
tified with the strongly correlated Auger feature of the 
CS theory. To underline this similarity further, we plot 
the average energy separation (A£^) and the intensity 
ratio (Ic/Im) between the correlated (satellite) and the 
uncorrelated (main) Auger spectral features as functions 
of U/W and /W'^, respectively in the insets (solid cir- 
cles) alongwith the results (open circles) from the CS 
theory. Clearly, both these show linear dependences for 
larger U. Between the two models we find that AE are 
similar even quantitatively. Ic/Im-, though qualitatively 
similar in both models, is larger in the multiband model. 

The spectra from partially-filled bands (Fig. 2) how- 
ever show qualitatively different behavior compared to the 
fully- filled limit. For brevity, we show the results only 
for some selected n and U/W values; we find that the 
results for other values of n and U/W are qualitatively 
similar. Figs. 2a and b show the results for U/W=0, 
1 and 2 in absence of any core-hole potential (Uc—O), 
while in Figs. 2c and d, we compare the results with 
and without Uc for the U/W=l case. In every case, 
the spectra appear a lot more complex than the simple 
and intuitive expectation of two groups of features aris- 
ing from uncorrelated and correlated final states in the 
Auger spectrum. Evidently, the spectra exhibit multiple 
groups of distinct and intense spectral features (marked 
1-5 in the figures for U/W—1 case) with as many as 
four spectral groups (2-5) appearing outside the energy 
region of the uncorrelated {U/W=0) spectrum (marked 
1). For example, the U/W=l case for n=3 (Fig. 2b) has 
features at about 1.5, 2.8, 4 and 5.3. Such results cannot 
be understood within the CS theory or in terms of the 
qualitative arguments presented in the introduction. Ad- 
ditionally, AE and Ic/Im exhibit no obvious dependence 
on U/W and U'^/W^ in contrast to the fully-filled case 
(Fig. 1). This is most apparent in Fig. 2a where the 
spectral feature 2 at about 1.3 moves to a lower energy 
with an increase of U/W from 1 to 2. In order to un- 
derstand the origin of such complex behaviors, we have 
analysed the character of final-state two-hole wavefunc- 
tions responsible for these various features. We illustrate 
schematically the dominant contributions to the wave- 



functions corresponding to each spectral group 1-5 in 
Fig. 2d in terms of electron and hole excitations with re- 
spect to the initial state. The screened, delocalized parti- 
cles (open and closed circles representing holes and elec- 
trons, respectively) are shown within the semi-elliptical 
bandwidth which is arbitrarily chosen to be half-filled 
for the purpose of clarity in the presentation. The local- 
ized hole states are shown in an atomic-like level. Thus, 
the spectral group marked 1 which corresponds to Auger 
transitions appearing within the energy interval of the 
uncorrelated {U/W=0) case corresponds to two screened 
holes in the band (schematic 1 in Fig. 2d). In each case 
the large intensity spectral group marked 2 was found to 
arise from an extra electron-hole excitation accompany- 
ing the generation of two screened holes (see schematics 
in Fig. 2d). The movement of the feature 2 to slightly 
lower energy with increasing U in Fig. 2a is due to the 
band-narrowing effects at larger U, as the process in- 
volves an excitation of an electron from the occupied to 
the unoccupied part. While the energy position of 2 is 
relatively insensitive to U (e.g. compare U/W^l and 2 
in Fig. 2b), the existence of this feature is entirely de- 
pendent on the presence of correlation in the initial and 
final states and consequently, this feature is completely 
absent for U /W=0. Feature 3 corresponds closely to the 
correlated feature discussed within the CS theory with 
two holes localized at the same atomic site, as shown 
by the schematic process 3 in the Fig. 2d. Thus, this 
feature moves to higher energy and rapidly gains inten- 
sity with increasing U (compare features 3 and 3' for 
U/W=l and 2, respectively in Fig. 2b). The very weak 
intensity feature 4 corresponds to an electron-hole exci- 
tation accompanying the generation of the two localized 
holes, whereas feature 5 at still higher energy arises from 
three holes localized at a site with an electron excited to 
a higher energy state. Though in these finite size calcu- 
lations, the features 2 and 4 appear as distinct peaks, 
in an extended metal such processes will have contin- 
uous energy spectra, thus substantially overlapping the 
corresponding main peaks, 1 and 3 respectively. These 
processes, however, will appear as distinct peaks in the 
case of insulators with finite band gaps. Obviously all 
these processes (2-5 with the exception of 3) are not 
describable within any theory for a fully-filled band, as 
it depends on multi-particle excitations possible only in 
partially- filled bands. 

The distinction between the partially-filled and fully- 
filled bands becomes all the more striking, when the in- 
fluence of the core-hole is included. It is evident that 
the core hole, existing only in the initial state, does not 
have any effect in the fuUy-fiUcd limit irrespective of the 
strength of the core- hole potential {Uc), since there can- 
not be any screening dynamics here. On the other hand, 
any reasonable value of Uc completely alters the Auger 
spectrum for any configuration away from the fully-filled 
limit, as we show in Figs. 2c and d illustrating typical 
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results for n=2 and 4 with and without the core-hole 
potential. In both these cases, correlation induced fea- 
tures 2 and 3 can be seen with substantial intensity for 
U/W=l in the calculated spectra in absence of a core- 
hole potential {Uc—0; thin solid line) with weaker features 
4 and 5 only for n— 4 case. For Ucj^ (thick solid line), 
however, the intensities in these features all but vanish 
completely, transfering almost all the spectral weight to 
the uncorrelated spectral range, 1. This observation of 
disappearance of the correlation-induced features in pres- 
ence of Uc is found to be robust for all reasonable values 
of U and Uc, as explicitly checked with extensive calcu- 
lations. Thus, one obtains here the paradoxical situa- 
tion where the pronounced manifestation of correlation 
effects within the valence band in terms of satellite fea- 
tures (2-5) is virtually wiped out due to the presence of 
another correlation effect, namely that between the core 
hole and the valence electrons. In order to understand 
these changes, we note that the spectra of the final state 
energy eigenvalues are identical for Uc—Q and ^ cases, 
since the core-hole exists only in the initial state. There- 
fore, the drastic modifications in the spectral features 
arise only from changes in the transition probabilities 
via the modification of the initial state wavefunction due 
to the screening response to the core-hole potential. The 
effect of the core-hole potential in the initial state is pri- 
marily to increase the local site occupancy by nearly 1 
electron to {n+1) as a screening response; analysis of the 
wavefunction character supports this view. Then, the 
subsequent Auger transition creates 2 holes which gener- 
ate primarily a {n-1) configuration at the local site, and 
therefore, suppress the correlation induced feature aris- 
ing from (n-2) local occupancy. It is to be noticed that 
there are some important changes in the line-shape of the 
spectral region 1 with correlation effects. In absence of 
the core-hole, there is a slight narrowing of this spectral 
region arising from band narrowing effects in presence of 
a finite U; with Uc ^ 0, there is a more significant narrow- 
ing in this spectral region (see Fig. 2d). This arises from 
the above-mentioned changes in the transition probabil- 
ities transferring weights preferentially to lower energies 
within the uncorrelated region. 

In order to ascertain the relevance of these results ob- 
tained from small, finite systems to strongly correlated 
extended solids, we have applied this method of calcu- 
lation to a representative case, namely the Auger spec- 
trum of a typical correlated oxide system, LaCoOs. The 
electronic structure of LaCoOs has been discussed ex- 
tensively in the past literature and there are reliable es- 
timates of the various interaction strengths already ob- 
tained from the analysis of several spectroscopic results 
within finite cluster many-body calculations involv- 
ing one Co and six oxygen atoms in an octahedral ge- 
ometry. In agreement with these estimates, we also use 
the same cluster with the charge-transfer energy A=2.0, 
pcicr=-1.8, A=A.2l, B=0.13, and C=0.64 amounting to a 



multiplet averaged Ueff of 4.5 eV. The calculated Auger 
spectrum which is the first of its kind is compared with 
experiment | [Tl| | in Fig. 3, indicating a remarkable agree- 
ment. Such an agreement between the calculated and 
the experimental spectra without the need to readjust 
parameter values firmly establish the relevance of the re- 
sults presented here. More importantly, it leads to a 
drastic reinterpretation of the experimental result. It has 
been believed that the most intense peak in the Auger 
spectrum in Fig. 3 at 8.8 eV for LaCoOs Q arises from 
Coulomb interaction driven localized two-hole final state, 
while the prominent shoulder at the lower energy (5 eV) 
arises from the screened and delocalized two-hole final 
states. An analysis of the wave-functions clearly shows 
that all the prominent features (i.e. at 5 and 8.8 eV) as 
well as the weak shoulder at about 13.3 eV arise from 
delocalized final states. We show calculated Auger spec- 
tra for LaCoOa with J7e//=4.5, 5.5 and 6.5 eV. Besides 
minor changes, the results are essentially identical inspite 
of the strong variation in U establishing our conclusions. 
Additionally, one can see the emergence of a broad and 
weak feature beyond 17 eV that moves to higher energy 
with increasing Uef / ; only such weak correlation induced 
satellite features survive in the final spectrum. In order 
to understand the origin of these spectral features, we 
first note that the Co d related single-particle density of 
states is distributed basically in two groups, arising from 
bonding-antibonding splitting of Co d - O p hybridized 
states [Q. Thus, the self-convolution of the single- 
particle density of states has three features corresponding 
to two holes occupying antibonding-antibonding (A-A), 
antibonding-bonding (A-B) and bonding-bonding (B-B) 
combinations at 4.7, 8.8 and 13.5 eV, establishing the 
origin of the experimental features at the same energies 
to be arising from A-A, A-B and B-B occupancies of the 
delocalized states. While the energetics can be explained 
by this procedure, interestingly the spectral shape cannot 
be described by a self-convolution of the DOS, due to the 
influence of the valence- valence and valence-core interac- 
tion effects modifying the transition probabilities. It is 
found that these correlation effects enhance the relative 
intensity of the peak at 8.8 eV. 

In conclusion, we have shown that the Auger spec- 
tra from strongly correlated, partially-filled systems are 
profoundly infiuenced by the presence of screening dy- 
namics and multielectron processes. Satellites in Auger 
spectra arising from correlation effects within the valence 
electrons is strongly suppressed by the presence of core- 
valence correlation effects leading to the surprising result 
of Auger spectra being dominated by uncorrelated spec- 
tral features inspite of the presence of strong interactions. 
This calculational method applied for the first time to a 
real system, LaCoOa shows a remarkable agreement with 
the experimental results, leading to a reinterpretation of 
Auger spectra from such strongly correlated systems. 
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I. FIGURE CAPTIONS 

Fig. 1 Calculated Auger spectra for U /W=0, 1 and 2 
for the fully-filled case. Insets show the variations in the 
relative energy position {AE) and the intensity {Ic/Im) 
of the correlated feature compared to the uncorrelated 
ones as a function of U and U'^, respectively for the multi- 
band (solid circles) and the single band (open circles) 
models. All energies are in the units of the bandwidth, 
W. 

Fig. 2 Calculated Auger spectra for (a) n=l, U/W=0, 
1, 2 and Uc=0; (b) n=3, U/W=0, 1, 2 and Uc=0; (c) 
n=2, U/W=0 and U/W=l with and without Uc^l.2 U; 
and (d) n=4, U/W=Q and U/W=l with and without 
C/c=l-2 U. Various groups of spectral features are marked 
as 1-5 and dominant contributions to the correspond- 
ing wavefunction are shown in panel d schematically. In 
panel b, high energy features for U/W=l are multiplied 
by 10 for clarity. All energies are in W . 

Fig. 3 Comparison of experimental and calculated Co 
'L^-VV Auger spectra from LaCoOs. The inset illustrates 
the changes in the calculated spectra with changing U , 
showing the existence of weak and broad correlation in- 
duced features at energies above ~ 17 eV. 
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